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ABSTRACT 

Using pulsewidth data for 872 isolated radio pulsars we test the hypothesis that pulsars evolve 
through a progressive narrowing of the emission cone combined with progressive alignment 
of the spin and magnetic axes. The new data provide strong evidence for the alignment over 
a time-scale of about 1 Myr with a log standard deviation of around 0.8 across the observed 
population. This time-scale is shorter than the time-scale of about 10 Myr found by previous 
authors, but the log standard deviation is larger. The results are inconsistent with models 
based on magnetic field decay alone or monotonic counter-alignment to orthogonal rotation. 
The best fits are obtained for a braking index parameter, n 7 « 2.3, consistent the mean of the 
six measured values, but based on a much larger sample of young pulsars. The least-squares 
fitted models are used to predict the mean inclination angle between the spin and magnetic 
axes as a function of log characteristic age. Comparing these predictions to existing estimates 
it is found that the model in which pulsars are born with a random angle of inclination gives 
the best fit to the data. Plots of the mean beaming fraction as a function of characteristic age 
are presented using the best-fitting model parameters. 
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1 INTRODUCTION 

One of the greatest of the many mysteries of radio pulsars is their 
evolutionary histories. As pulsars age, powerful electromagnetic 
torques act to increase the rotation period, P, of the strongly mag- 
netized neutron star, causing the spin-frequency, v = 1/P, to de- 
crease over time, a phenomenon known as spin-down. This can be 
described by the equation 

v = -Kv n , (1) 

where n is a parameter known as the braking index, with n — 3 
for a magnetic dipole rotating in a vacuum. The spin-down torque 
acting on the star is proportional to v, and it decays over time as v 
decreases. 

There is also evidence that K is not constant, but rather 
is also changing over time. The surface dipole magnetic field 
strength at the magnetic equator, assuming that the spin and mag- 
netic axes are orthogonal, is conventionally given by B SUI f — 
3.2 x 1 19 (PP) 1/2 G; at th e magnetic poles the field strength is 
2B surf dLvne & SmitifeoOfj) . P aur f tends to be bigger for younger 
pulsars than for older ones having a large spin-down age, t s = 
P/ [(n — 1)P] . This in turn suggests that K tends to reduce with 
increasing age, and much conjecture remains about the origin(s) of 
this evolution. 

Many authors have argued that the main contributor to 
the i\"-evolution is the progressive alignment of the spin and 
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magnetic axes: e.g. |j_ :ancly 
dl986b, iLvne & Manchester! 
i Kuz'min & WvJ j 1992b. ICandvl 
| Tauris & Manchester! dl998h and 



1988) 
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Pandev & Prasa< 



1991) 



1996) 



Tauris & Manchester (1998) and Weltevrede & Johnston (2008) 
But iMcKinnod jl993h and iGil & Had Jl996h have "argued 
favour of a random distribution of the inclination angles between 
the angular velocity and magnetic axes and against magnetic 
alignment. Others have favoured magnetic field decay (e.g. 
iNaravan & Ostriker|[l990r) . while there are also authors who argue 
for continuing counter-alignment to eventual orthogonality of the 
axes, via an electromagnetic torque exert ed by magnetosph eric 
electric currents fl owing on ope n field lines dBeskin et al Jl 1984b . 

According to Jonesj dl976b , pulsars initially counter-align (via 
a strongly temperature-dependent dissipative torque in the fluid in- 
terior) to reach the orthogonal rotator state after about 10 3 yr, where 
they remain for 10 4 -10 5 yr (while the dissipative torque decays 
as the interior cools), after which alignment (via electromagnetic 
torque) occurs. 

If K were constant, the braking index n would be equal to the 
apparent braking index, 



vv jv 



pp/p 2 



If v can be determined from observations of a pulsar, then n a 
be computed. Howev er, if K is varying, then in general 



(2) 
, can 

in general n app ^ 

n. So me authors (eg. Johnston & Galloway 1999'; lTauris & Konar] 
l200lh have argued that, for pulsars of moderate age, n app exceeds 
the value of 3 that corresponds to magnetic dipole radiation. This 
could be taken to imply that either magnetic alignment or field de- 
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cay must be occurring. However, it is likely that recovery from un- 
seen pulsar glitches is the pr incipal contribution to measured vari- 
ations in P (Wa ns et alj200lh for these pulsars. 

In this paper, we shall use the p ulsar evolution model s from 
ICandv & Blaiifd 1983L hereafter CB83). ICandv & Blain ( fl986l here- 
after CB86) and lCandvl d 19931 , hereafter C93), which all invoke a 
progressive alignment of the spin and magnetic axes. The CB83 
model develops the pulsewidth-age relation on the assumption that 
all pulsars move along the same evolutionary track, differing only 
in age and orientation with respect to Earth. The CB86 model re- 
laxes this assumption and allows a distribution of alignment time- 
scales, as well as a time-varying relationship between characteristic 
and actual age. The C93 models additionally allow a distribution of 
initial inclination angles. All of these models predicted a minimum 
in the mean pulsewidth as a function of characteristic age and pro- 
vided a good fi t to the available data, consis ting of 293 pulsewidth 
values from the lManchester & Tavloti 1 1981) catalogue. 

Here we use t he ATNF (Australia Tel escope National Facility) 
pulsar catalogu fl( Manch ester et al.l2 005). version 1.35, to investi- 
gate pulsar evolution through magnetic alignment using the Candy- 
Blair models. This contemporary catalogue contains much larger 
data sets of 872 and 1420 isolated radio pulsars with 10 per cent and 
50 per cent intensity pulsewidth values, Wio and W50. It also con- 
tains many more old pulsars than was previously the case, which 
tend to have lower radio luminosities. It thus allows a more pre- 
cise analysis of pulsar evolution. As a consistency check, we also 
examine separately those pulsars with pu lsewidths measured dur- 
ing the Parkes m ultibeam pulsar survey dManchester et alj|200ll ; 
lLorimer et al.l2006h . giving Wio and W 5Q data sets of 377 and 934 
pulsars respectively. 

The paper is structured as follows. Section[2]presents a sum- 
mary of the theory of alignment between the spin and magnetic 
axes; Section|3]gives a mathematical description of pulsewidth evo- 
lution; in Section|4]we least-squares fit the evolutionary pulsewidth 
models to the Wio and W50 data as functions of characteristic age; 
in Section[5]we further analyse these results, in particular co mpar- 
ing them to t he inclination angle estimates of lRankinl dl993bl) and 
lGoulJdl994h ; and conclusions are presented in Section|6] 



2 MAGNETIC ALIGNMENT THEORY 
2.1 Evolution in true age 

The Candy-Blair pulsar evolution models are based on two effects. 
One is progressive alignment of the spin and magnetic ax es, for 
which the formula describing the alignment phase of the Ijonesl 
d 19761) model is adopted: this has an exponential decay of the sine 
of the magnetic-spin inclination angle, a, from its initial value etc,: 



p(t)<xP->(t) : 



(4) 



sina(t) = exp(— t/r) sinao . 



(3) 



where t is the pulsar's age and r is the alignment time-scale. The 
alignment has a simple electromechanical origin: the electromag- 
netic radiation emitted by an oblique rotating dipole results in a 
torque on the star that causes the angular velocity axis to migrate 
through the neutron star toward alignment with the magnetic axis. 

The second effect is progressive narrowing of the emission 
cone, described by a power-law dependence of its half-angle, p, on 
the rotation period, P: 
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with 7 a positive constant having a value between 1/3 
and 2/3 dGunn & Ostriketlll970l: :l Ruder man & Sutherland [1971: 



iLvne & Manchestedll988l ; lRankinlll993d) . 



For a pulsar with constant rotational inertia and fixed magnetic 
moment vector, the rotation period ev olves according to pp n ~ 2 = 
constant. It was pointed out long ago (Phinnev & Blandfordi ll98lh 
that this simple rule is incompatible with the observed distribu- 
tion of pulsars in the P— P plane, conflicting with stationary 'flow' 
through that plane to pulsar 'death' at advanced ages. However, 
alignment produces an effective reduction in the magnetic moment, 
reducing the torque in proportion to sin 2 a; so, in the Jones m odel, 
the period evolves according to dPhinnev & Blandforj|l98lh 



PP n ~ 2 oc exp(-2t/r) , 



with n (the braking index) constant. 
Equation ((5) integrates to give: 



+ (Po r ' 



)exp(-2t/T) 



(5) 



(6) 



where Po is the pulsar's initial period and Poo is its limiting pe- 
riod as t — > 00. Assuming that pulsars eventually slow substan- 
tially from their initial spin period, Po <C Poo, we can approximate 
equation ((6) by 

P(t) « Poo[l - exp(-2i/r) f'^ . (7) 

Using for P(t) in $4$ for p(t) yields a simple approximation for 
the evolution of p, expressed directly in terms of t: 



p(t)^Poo[l-e X p(-2t/r)]^ /( "- 1) . 



(8) 



where poo is the limiting emission cone half-angle as t — > 00. 

We note that, because Po has effectively been taken as zero in 
going from to <[7j for P(t), equation ([Sj( is not applicable to the 
very youngest pulsars: for t <C r/2 it reduces to 



(9) 



which diverges as t — > 0. For the typical parameter values (see 
Section 1431 Table [3l 7 = 1/2, n = 2.3, r = 1 x 10 6 yr, and 
poo = 2 — 5°, equation l[9} gives p(t) ~ llpoo ~ 22 — 55° fort ss 
1000 yr. As only one known pulsar is this young, the divergence is 
outside the age range we can analyse, and the approximation d8} 
for p{t) is valid for our purposes. 



2.2 Evolution in characteristic age 

A pulsar's spin-down age, t 3 = P/ Un — 1)P] , depends on n, 
which is going to be a variable in our fitting procedure below. In- 
stead of the spin-down age, from here on we use the the character- 
istic age which is defined to be 



tc 



PI (2P) . 



(10) 



The characteristic age is independent of n and corresponds to the 
spin-down age in the case n = 3. Its use enables us to show fittings 
for different n on one plot using the same characteristic-age data. 
Inserting equation for P(t), we find that t c is related to the 
pulsar's actual age, t, by: 

t c (t)/r = [(n-l)/4][exp(2Vr)-l] (11) 
or exp[2t(t c )/r] = 1 + (4t c )/[(n - 1)t] . (12) 

Note that factors of 2/(n — 1) effectively rescale t c to t s in these 
equations. 
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It follows from equation dl lb that t c ~ (n — l)t/2 for rel- 
atively young pulsars (young compared with the spin-down time- 
scale); for old pulsars, t c /r increases as an exponential function 
of t/r. The characteristic and actual ages coincide for relatively 
young pulsars in the dipole rotator case, for which n — 3. 

The t c -t relation enables equations l[3j, Q and l[8j for a, P 
and p to be re-expressed in terms of t c - which has the advantage 
of being directly determined from the measured quantities P and 
P - instead of in terms of t: 

sina(i c ) = (sina ) [1 + 4t c /(n - l)rp 1/2 , (13) 
P{t c ) = Poo [1 + (rc - l^/^n 1 /'™- 1 ' , (14) 
and p(t B ) = Poo [1 + (n - l)r/4t c ]^ /( "- 1) . (15) 

Equation J 1 5t describes an emission-cone half-width that de- 
creases quite rapidly with time for young pulsars; consequently, the 
observed pulsewidths should decrease with age for the youngest 
pulsars. The alignment of the spin and magnetic axes, described by 
equation d!3t , should cause an increase in the observed pulsewidths 
for older pulsars, as the pulse will occupy a larger fraction of a com- 
plete rotation. The combination of these two processes should lead 
to a minimum in the observed angular pulsewidths for pulsars of 
moderate ages (CB83). 

Note that taking the derivative of equation $6$ and using equa- 
tion l |10t for t c gives (cf. equation l|2}) 



and equation d!6t gives 



= n + U c /t , 



(16) 



showing that the apparent braking index increases with character- 
istic age. 

2.3 The parameters of the theory 

Equations d!3t . d!4t and d!5t - expressing the evolution of the 
magnetic-spin inclination angle, the pulsar period and the emis- 
sion cone half-angle in terms of a pulsar's characteristic age - can 
conveniently be re-written in terms of a normalized dimensionless 
characteristic age: 

T c = 2t c /[(n - 1)t] , (17) 

which reduces to t c /r in the dipole rotator (n = 3) case. Thus: 

sina(T c ) = (sina )(l + 2T c r 1/2 , (18) 

P{T C ) = P 00 [l + (2T c )- 1 r 1/(n ' 1) , (19) 
and p(T c ) = p 00 [l + (2T c )- 1 ]^ /(n - 1) . (20) 

In the equations used for the pulsewidth modelling below, 
namely J 1 3 b and d!5l l, or d!8b and d20t . for a(t c ) and p(t c ), n ap- 
pears only in the two combinations 7/(71 — 1) and (n — l)r. Hence 
there are only four basic independent parameters in the Candy-Blair 
alignment models. Later we will s et 7 = 1/2, consistent with the 7 
values found by other authors (e. g. lRankinl 1 9 93a ; McKinno nl 19931 ; 
Gould i 19941) . Hence, it is convenient in this study to define the four 
independent parameters to be Qo, Poo, n-y and t 7 , where 

n 7 = 1 + (n - 1) /(27) , and (21) 
r 7 = 27 r, (22) 

since these last two parameters are normalized to the 7 = 1/2 case. 
In terms of these parameters, the two combinations 7/ (n — 1) and 
(n — 1)t become 



7/(71-1) 

(n- 1)t 



1/ [2(n 7 - 1)] , and 
(ti 7 — l)r 7 , 



(23) 
(24) 



1 = 27 In. 



1 + — 

r 7 



(25) 



3 PULSEWIDTH EVOLUTION DESCRIPTION 
3.1 Mean pulsewidths 

For an individual pulsar, the observed angular pulsewidth is given 
as a function of a, p and ( by 



W(a,p, — 2arccos 



cos p — cos a cos £ 
sin a sin £ 



(26) 



where £ is the angle between the observer 's direction and the pul- 
sar's spin axis dManchester & Tavlojl 19771 . p. 218). The strong de- 
pendence on the orientation of the observer with respect to the spin 
axis leads to a wide scatter in the pulsewidth values, so we need to 
take a mean pulsewidth for pulsars with certain values of a and p. 

The mean pulsewidth is obtained by integrating d26t over the 
angular extent of the emission beam: 



(W){a,p) = 



W(a,pX)P(0d(. 



(27) 



a-p\ 



in which P(£)d£ is the probability that the angle between the spin 
axis and the observer's direction is in the range £ to £ + d(\ The 
integration is over all angles £ for which emission is both directed 
toward the observer and seen as pulsed. 

The limits of integration in d27t for {W)(a, p) are explained 
by means of Fig. Q] For a pulsar that is not too close to align- 
ment, such that the emission beam does not contain the rotation 
axis (Fig.Q] left), a > p and the beam is intercepted for a part (but 
not all) of the rotation period if £ is in the range a — p to a+p. For a 
pulsar that is close to alignment, with the emission beam containing 
the rotation axis (Fig.Q] right), a < p and the range £ < p — a is 
excluded, as some part of the beam is then always directed toward 
the observer, assuming the emission cone to be filled. 

We take a ^ 90°; the choice of first or second quadrant for 
a is a matter of convention, switching it between hemispheres be- 
ing equivalent to reversing the sense of the neutron star's rotation. 
The beam, however, can extend into the other hemisphere, which it 
will do if a + p > n/2, as expected for young pulsars with wide 
beams. If a pulsar produces observable radio emission from both 
poles, then this will be observed as a main pulse and interpulse, 
provided that the pulsar is nearly orthogonal or the pulsar beams 
are wide. However, the mean pulsewidths referred to here are for 
emission from only one of these poles, not the combined main pulse 
- interpulse width. 

On taking pulsar spin axes to be randomly directed with re- 
spect to observers, and normalizing over the angular extent of the 
emission cone, the probability distribution for £ becomes: 



P(C)dC=(sinCdC/2)// fj (a,p) 



(28) 



where /(, is the 'beaming fraction', which measures the fraction of 
the sky swept out by the beam. Equation d28t is the ratio of the 
surface area of a ring on a sphere, corresponding to observers in 
the range £ to £ + d£, to the surface area of the sphere swept out 
by the beam. 

The surface area of the part of a sphere between co-latitudes 
61 and 62, normalized to the surface area of a sphere, is 



sin [(0! + 6> 2 )/2] sin [(0 2 - 6>i)/2] 



(29) 
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Figure 1. Pulsar emission relative to a distant observer, assuming a circular 
cone of emitted radiation, with the magnetic-spin inclination angle a either 
greater (left) or less (right) than the emission cone half-angle p. In the latter 
case, the emission cone encloses the spin axis, and there may be little or no 
modulation; in particular, to observers in the cone £ < p — a, the emission 
cone is continuously visible. The left-hand diagram depicts the usual case 
of a pulsar with moderate values of a and p, typical of mature-aged pulsars. 
The right-hand diagram depicts a pulsar with a thin beam and nearly aligned 
spin and magnetic axes, typical of old pulsars in the Candy-Blair theory and 
the Jones model. 



Inserting 9i = \a — p\ and 62 = a + p gives the beaming fraction: 

fb(a, p) = sin a sin p , (30) 

applicable to both situations depicted in Fig.Q] 

We note that lTauris & Manchester! ( 119981 , eq. (7)) use a differ- 
ent beaming fraction for the case a < p, namely 1 — cos(a + p), 
which results from integration over to a + p. As they remark, 
there may be little or no modulation when a < p, and we prefer to 
exclude the 'continuously seen' cone, £ < p — a (Fig.[T| right). 



3.2 Evolution modelling 

The evolutionary time dependence of a and p, equations ([3} and 
l[8}, results in a mean pulsewidth, (W)(a,p), that is a function 
of time only, (W)(t). This function is to be treated as a density, 
(W)(t) dt giving the mean pulsewidth of pulsars aged between i 
and t + dt. The age of the pulsars will here be measured in terms of 
their log characteristic age, log t c , and will be binned with respect 
to that variable. Under the change of variables t — > i(logi c ), we 
have 



(W)(t)dt = (W)(t(\ogt c )) 



dt 



d (logi t 
(W) (log t c )f t d (logt c ) 



-d (logi c 



The factor ft (i c ) is given by 
di 



d log i c 



rln 10 /, (n- 1)t 
1 + A + 

4ic 



(31) 
(32) 

(33) 



where the t-t c relation, equation d!2l >. has been used. So f t (i c ) in- 
creases monotonically from (In 10)2i c / (n — 1) for young pulsars, 
flattening to (In 10)r/2 for old ones. 

To allow for variations between the evolutionary histories of 
different pulsars, CB86 and C93 relaxed the assumptions that all 
pulsars are born with the same inclination angle and have the same 



alignment time-scale, an and r respectively. Instead, these are re- 
placed by a probability distribution function, P(ao), for initial in- 
clination angles (see Section |33l , and a log-normal distribution of 
alignment time-scales: 



P(logr) = 



27T<7 lo . 



■ exp 



°"log 



(34) 



where pi og and rji og are the mean and standard deviation of the 
logarithm of the time-scale distribution, logr. If the normalized 
parameters n 7 and r 7 that were introduced in equations d2 1 1 — d241> 
are used, then r 7 has a log-normal distribution with mean 



= Mlog + log (27) 



(35) 



and standard deviation ai og unchanged. 

The mean observed pulsewidth as a function of characteristic 
age is now given by 



(W) (logic) = 



00 vr/2 



{W)(a,p)f b ftP(a )P(logr)da 



dr 





in which the normalization factor, A, is 

rn/2 



A(logt c ) = 



(i 



dr 



f b ftP(a a )P{log T )da — 



(36) 



(37) 
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The beaming fraction, f b , must be included as a weighting factor 
in equation l |36t to take into account the effect of beaming on the 
probability of observation, and those weights are normalized by 
the factor A^ 1 (the mean pulsewidth is only computed for those 
pulsars that are observed). The aim here is to fit equation l |36t for 
(W)(logi c ) to the observed data. 
In Appendix|A]it is shown that 



{W)(a,p)f b = 



[cos(p - 



- COS p\ 7T 
a) — COS p] 7T 



for 
for 



a ^ p , 
a < p . 



(38) 



Inserting this result into equation l |36t for (W) (logic) allows a 
considerable improvement in the numerical performance of the 
curve-fitting. 

At a given characteristic age, Equation l |15t implies that 

p(t c ) > 90° if a pulsar were to have r > r cut , where 



Tout 



(tc) 



4ic 



4i e 
n- 1 



90° 

poo 



(n-l)/j 



- 1 



, for poo < 90° 



90° 

Pc 



(39) 



(40) 



This unphysical possibility arises from the approximation made 
earlier that Pn / Poo — * 0. In reality a pulsar with r > r cut must 
be born with a characteristic age greater than the i c used in equa- 
tion l |39t . Therefore integrals over r are terminated at an upper limit 
of 0.99r cut , where the factor 0.99 is included to assist numerical 
convergence. 



3.3 Three initial alignment models 

We use three different initial alignment models, P(cto), from C93 
that involve different distributions of the initial inclination angles. 
Model I assumes that all pulsars start with their spin and magnetic 
axes mutually perpendicular: an = 90°. This corresponds to the 
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Ijonesl model following the initial brief phase of counter-alignment 
to orthogonality. Model II uses the probability distribution 

P(a ) = (2 - V2)" 1 sin(a /2) , (41) 

which assumes pulsars are born with a random distribution of in- 
clination angles from to 90° . Model III uses 

P(ao) = (4/tt) cos 2 qo , (42) 

which assumes pulsars tend to be born with inclination angles 
closer to than would occur by chance. At the outset there is no 
particular reason to think that Model III might be approximately 
correct; it is included here for contrast with Models I and II. 



4 FITTING TO DATA 
4.1 Data selection 

Both 10 per cent and 50 per cent intensity pulsewidth data, Wio 
and W50, are used in this analysis of pulsewidth e volution. The data 
were retrieved from the ATNF pulsar catalogue ( Manchester et al.l 
2005), which includes results from several surveys at around 400 
and 1400 MHz. For this study, all the binary, millisecond, anoma- 
lous X-ray, rotating radio transient (RRAT) and globular cluster 
pulsars have been removed from the data set, because these pul- 
sars are considered to follow different evolutionary histories from 
those of the more common isolated radio pulsars. The condition for 
excluding millisecond (or low magnetic field) pulsars was chosen 
to be B SUI{ sC 4 x 10 9 G. 

Pulsars with interpulses whose Wio or W50 values potentially 
include both the main pulse and interpulse were also removed. To 
do this the complete census of the interpulse puls ar population 
recently published by Welt evrede & Johnston! d2008h was used. It 
identified 27 pulsars with interpulses, and our Fig. [2] shows his- 
tograms of the existing Wio and W50 values for them: it is evident 
that the distributions can be treated as bimodal with one maximum 
well below 90° and the other at around 180°. For the purposes 
of this study it was assumed that a value of Wio or W50 > 140° 
includes emission from both poles of the pulsar (that is both the 
main pulse and the interpulse), and so should not be included in the 
study. This resulted in 6 pulsars being removed from the Wio data 
set leaving a net total of 872 pulsars; and 3 interpulse pulsars being 
removed from the W50 data set, leaving a net total of 1420 pulsars. 
These two data sets are referred to hereafter as ATNF Cat Wio' 
and ATNF Cat Wso' ■ 

In a mean pulse profile the spacing between components usu- 
ally increases, and with it the whole profile width, as the radio 
frequency, /, decreases. For ins tance, by mea s uring the profile 
widths of 6 pulsars up to 32 GHz. lXilouris et al] ( 1996) found that 
Wso = bo + bif~ p where p, bo and bi are constant for each pul- 
sar, with 0.3 < p < 0.9. The Wio and W50 pulsewidths tabulated 
in the ATNF catalogue are not all measured at the same observ- 
ing frequency. The frequency dependence of the pulsewidth there- 
fore introduces some statistical noise into the ATNF Cat Wio' and 
ATNF Cat W50' data sets. As a consistency check, we have also 
examined the subsets consisting of only those pulsars whose cata- 
logued pulsewidth val ues were measured in the Parkes Multibeam 
Survey at 1374 MHz ^Manchester et al.l20oil ; lLorimer et al.l2006t) . 
These subsets are referred to hereafter as 'Parkes MB Wio' and 
'Parkes MB W50', and contain 377 and 934 pulsars respectively. 

Both the Wio and W50 characteristic age-pulsewidth data dis- 
tributions naturally lend themselves to binning in log t c , and it was 




50 100 150 200 250 



Pulse Width (deg) 

Figure 2. Frequency histograms of the H^io (b lack) and W50 (grey) 
pulsew idths for the interpulse pulsars listed by Weltevre de & Johnston! 
12008); the bars are stacked for clarity. It is evident that the data distri- 
butions can be treated as bimodal. Widths greater than 140° were assumed 
to have included the interpulses and were excluded from the Wio an d W50 
pulsewidth data sets used in this paper. 

also found numerically optimal to fit the mean pulsewidth, equa- 
tion J36t . as a function of log t c rather than t c , i.e. as ( W) (log t c ). 
The Wio or W50 data set is then binned to give N data points, 
{ (log t c ) . , Wi}, i = 1, . . . , N. It is desirable to make the bins as 
narrow as possible in order to give maximum resolution in logt c . 
This is because if a bin i extends over a large range of log t c in a 
region where the underlying function, here (W)(logt c ), has a sig- 
nificant slope, then a part of the uncertainty in the estimated value 
of the mean pulsewidth, Wi, will result from the change in the 
underlying (W)(log t c ) across the bin. On the other hand, it is de- 
sirable to include as many points as possible in a bin to reduce the 
uncertainty in the Wi estimate of (W) [(log i c ) .] , and also to en- 
sure that the errors in the estimates are approximately normally dis- 
tributed so that standard data-analysis procedures can be applied. 
To attempt to balance these competing effects we binned the data 
so that each bin contained about 30 pulsars - the minimum num- 
ber of points to ensure that the errors are approximately normally 
distributed. 

Least-squares fittings are made to both the Wio and W50 data 
sets. Even though there are more data points in the W50 than the 
Wio data sets, the Wio is superior for our pulsewidth fitting be- 
cause the W50 set suffers more from confusion. This is because 
a W50 measurement will include outlying components in a profile 
only if they are above 50 per cent of the maximum intensity. A Wio 
measurement is more likely to include the outlying components and 
so will give a more consistent estimate of the total pulsewidth. A 
histogram of the ratio of W10/W50 is found to be unimodal, with 
a peak at about 2. However, it extends to ratios above 10, demon- 
strating that the Wio in those cases is detecting outlying emission 
missed by the W50 measurement. The trade-off in relying on Wio 
values is that the Wio measurement is more difficult to make than 
the W50, and so there are fewer pulsars with existing Wio measure- 
ments. 

The data sets considered in this paper are summarized in Ta- 
ble Q] For each data set it gives the name of the set and the num- 
ber of pulsars it contains, and describes the binning of the data. 
Fig-H]shows the ATNF Cat data sets and compares them to the cor- 
responding Parkes MB data sets. We see that the ATNF Cat and 
Parkes MB data sets have very similar distributions with respect to 
characteristic age. 
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Table 1. Summary of the four data sets considered in this paper. 



Dataset 


Size 


Binning 


Parkes MB Wio 


377 


bins 1-13: 29 pulsars/bin 


Parkes MB W50 


934 


bins 1-14, 19-31: 30 pulsars/bin 






bins 15-18: 31 pulsars/bin 


ATNF Cat Wio 


872 


bins 1-14, 17-29: 30 pulsars/bin 






bins 15-16: 31 pulsars/bin 


ANTF Cat W 50 


1420 


bins 1-19, 30^17: 30 pulsars/bin 






bins 20-29: 31 pulsars/bin 




10 r 

4 5 6 7 8 9 

50 I 



40 - 




log(f c /years) 

Figure 3. Comparison of the four mean-pulsewidth data sets considered in 
this paper, as summarized in Table[T] The Wio and W50 'ATNF Cat' data 
sets are shown in black (connected by a dashed black line), while the two 
corresponding 'Parkes MB' data sets are shown in grey (connected by a 
solid grey line). The error bars are \a standard errors for the binned data. 



4.2 Fitting procedure 

Each model of initial alignment has four free parameters: p^, n 7 , 
p 7 ,i og and criog. The aim is to determine these four parameters 
by least-squares fitting the model (l4 y )(logt c )-curve to the binned 
data. The mean pulsewidths, Wu computed for each bin i, have 
normally-distributed errors, and so the uncertainties for the fitted 
parameters can be found by com puting the parameter covariance 
matrix, Ski (e.g. |Pres7 et al. 1992, Ch. 15). This matrix can in turn 
be used to compute a parameter correlation matrix, rtu for the pa- 
rameter estimates. 

Proceeding in this way it was found that least-squares esti- 
mates of the four parameters had large uncertainties, and hence 
were of questionable statistical merit. These large uncertainties 
were found to result mostly from a high correlation between the 
four parameters, with the off-diagonal values of the parameter cor- 



G. Blair 

ATNF Cat W w Data 

23 r 



2.2 - 




1.5 2.0 2.5 3.0 3.5 4.0 4.5 



n 7 

Figure 4. Plot of Xrcd as a f unctlorl of n 7 , for each of the three models of 
initial alignment fitted to the ATNF Cat Wio data: Model I (solid black), 
Model II (dashed), and Model III (dotted). Each model exhibits a minimum 
in Xrod at n f aroun( l 2.3. 

relation matrix all being around ±0.98. Furthermore, very different 
sets of parameter values could give very similar fittings to the data. 

To tackle these problems we fixed n 7 at specific trial values, 
leaving three free parameters: p^, p 7 ,i og! and ai og . The cross- 
correlation between the parameters was significantly reduced. Cor- 
respondingly, the uncertainties in the parameters were also signif- 
icantly reduced and, given the assumptions, provided statistically 
significant constraints. The reduced \ 2 (Xrcd = X 2 /dof) for these 
three-parameter fits were then examined as a function of the pre- 
set parameter n 7 ; a minimum in this function would indicate that 
a particular value of n 7 was giving a superior fit. In that case a 
least-squares fit was computed by allowing all four parameters (in- 
cluding n 7 ) to vary, and the solution was checked against the three- 
parameter fit for consistency. 

We recall that the parameters n 7 and p 7 ,i og are the values 
of n and /zi og assuming that 7 = 1/2. More generally, for fitted 
values of n 7 and pi 7 ,i og , the relationships between different values 
7, n and /ii og are found via equations l !2 1 1 and J351 > to be n = 
27 (n 7 — 1) + 1, pi og = p 7 ,i og — log (27), and hence 

n — 1 
n 7 — 1 

4.3 Fits to the data 

The three-parameter fits showed clear minima in x? e d ( n j) at n i ~ 
2.3 for all of the datasets, excluding the ATNF Cat W50 data. For 
instance, Fig.|4]plots Xied i n i) f° r Models I, II and III fitted to the 
ATNF Cat Wio data. 

Fig-HJshows three-parameter fitted curves of (W)(logt c ) for 
Model II with n 7 set at values between 1.5 and 4; equivalent plots 
for Models I and III are similar. In each case setting n 7 to around 
2.3 and least-squares fitting the other three parameters gives the 
best fit to the data, in particular at low characteristic age. Perform- 
ing the same procedure with n 7 higher or lower than this value 
results in a curve that goes to lower pulsewidth than the measured 
value at low log t c . The figure shows for the ATNF Cat VK50 data set 
that even though x? c d does not have a clear minimum at n 7 w 2.3, 
the n~f ~ 2.3 fit appears to best explain the low-t c mean 147 50 
pulsewidths. 

Table [2] gives the measured values of n app for six pulsars (ta- 
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ATNF Cat W 10 Data & Model II Parkes MB W 10 Data & Model II 




log(? c /years) log(r ( ./years) 



Figure 5. Mean pulsewidth, (W), as a function of log characteristic age. The data sets are given in the title of each figure; they are the same as those plotted 
in Fig.[3]and are summarized in Table[T] Error bars are lcr standard errors for the binned data. Mean angular pulsewidths for Model II are least-squares fitted 
to the binned data for each of the following values of n 7 : n 7 = 1.5 (dashed), n 7 = 2.3 (solid), n 7 = 3.0 (dot-dashed), and n 7 = 4 (dotted). The vertical 
scales are different for the Wio and Wbo data sets. Note that for the Parkes MB W50 data set, the Model II curves did not converge for n 7 = 1.5 or n 7 4. 
In that case, the dashed curve and the dotted curve represent n 7 = 1.7 and 3.95 respectively. In each plot, except the ATNF Cat W50 plot, the solid n 7 = 2.3 
curve is approximately the best-fitting curve to the data across the different values of n 7 , having the lowest value of Xred • None of the curves in the ATNF 
Cat W50 plot was clearly distinguished by a minimum in X? e d ' 



ble 2, p. 1291. lLivingstone et alj2 006. and references therein). The 
uncertainties in the last digit of n app are shown, being la confi- 
dence intervals. The characteristic ages are also listed, indicating 
that these pulsars are all very young. Except for the Vela pulsar, the 
determinations of n app were all made by standard timing analyse s 
that obtained phase coherent solutions dLivingstone et alj l2007h . 
This approach could not be used for the Vela pulsar because of its 
large glitches, so v was determined from data 150 days after each 
glitch and then extrapolated back to the glitch epoch; v could th en 
be determined from the change in v over time jLvneet all 1 996). 

It is noteworthy that the mean value of n apP in Table [2] is 
2.42 ± 0.05, compared to the values of n 7 ~ 2.3 that give the 
best fits of the Candy-Blair models to the pulsewidth data. On the 
grounds of equation ((25) we expect n app w n 7 for 7 w 1/2 and 
t c <C t 7 . Fig. [5] shows that the constraint of n 7 « 2.3 is based 
primarily on young pulsars, but a much larger sample than the six 
given in Table [2] 

Table [3] gives the parameter values obtained from the three- 
parameter fits to the four data sets. The error ranges are l<r confi- 
dence intervals, and the degrees of freedom (dof) for each data set 
are given. For each fitting the parameter n 7 is fixed at the value at 



Table 2. Measured values of the apparent braking index, n : 



Pulsar 


tc (yr) 


f-app 


Reference 


J1846-0258 


723 


2.65(1) 


Livingstone et al. (2006) 


B0531+21, Crab 


1240 


2.51(1) 


Lvne et al. (1993) 


B 1509-58 


1550 


2.839(3) 


Livingstone et al. (2005) 


Jl 119-6127 


1610 


2.91(5) 


Camilo et al. (2000) 


B0540-69, LMC 


1670 


2.140(9) 


Livingstone et al. (2005) 


B0833-45, Vela 


11300 


1.4(2) 


Lvne et al. (1996) 



which Xvcd 1S minimum, except in the case of the ATNF Cat W50 
data where no clear minima were found, and so n 1 was fixed at 
2.3 for comparison with the other fits. For each data set, all three 
models predict similar /U 7 ,i g and a\ og values. The only signifi- 
cant parameter difference between the three alignment models is 
the predicted limiting emission cone half-width, poo. Fig. [6] shows 
curves of (VK)(log t c ) for Models I, II and III using the parameter 
values listed inTable[3] It is evident that the pulsewidth data alone 
can not be used to distinguish between the three models of initial 
alignment. 
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Figure 6. Mean pulsewidth, (W), as a function of log characteristic age using the fitting parameters given in Table|3] Mean angular pulsewidths for Models I 
(dashed line), II (solid) and III (dotted) are least-squares fitted to the binned data. 



The figures show that following a rapid decline in pulsewidth 
for young pulsars, there is a minimum in the pulsewidth for pul- 
sars of moderate age followed by an increase for older pulsars: this 
is a clear signature of spin-magnetic alignment (CB83). The min- 
imum occurs at t c ~ 10 6 ' 9 yr in all three models for all data sets, 
except the Parkes MB Wio data set where it occurs slightly earlier. 
The mean alignment time-scales obtained from the fitting process 
are around 10 6 yr which is less than previous estimates of around 
10 7 yr. The log deviations of the time-scales range from 0.71 to 
1.21, and are all larger than the value of 0.25 found by both CB86 
and C93, which suggests that the more sensitive recent surveys have 
detected pulsars with a broader range of parameters. 

As discussed near the end of Section |4~T1 Wio is superior to 
Wso for pulsewidth fitting. Hence in the following we restrict our 
attention to the Ww fits, in particular the ATNF Cat Ww fit that 
was generated from a larger data set. 



5 ANALYSIS AND IMPLICATIONS 
5.1 Magnetic field decay model 

As an alternative to magnetic alignment, many authors have pro- 
posed pulsar evolution through magnetic field decay. This will have 
a different effect on pulsewidth evolution. Suppose that the mag- 
netic field decays exponentially, according to 



B = S exp(-t/r d ) , (44) 

where Td is the field decay time-scale, and that a — ao remains 
constant (no alignment). Equation $5$ for P(t) still applies, so we 
still have equations {8]l and J 1 5 b for p(t) and p(t c ). 

Fig. |7] shows curves of (W)(\ogt c ) for a constant a = Qo 
(using Models I, II and III for the ao distribution) to illustrate the 
effect of magnetic field decay on pulsewidth evolution. The fig- 
ure shows the same ATNF Cat Ww binned data as before and the 
curves have been fitted assuming that log Td has a normal distribu- 
tion with mean fid and standard deviation ad. It was found neces- 
sary to fix ad in order to get a convergent fit, and the small value 
of ad = 0.1 was required in order to increase the size of fid to 
realistic values. 

Qualitatively, the effect of field decay alone on pulsewidth 
evolution is a progressive decrease in pulsewidth over time, as seen 
in Fig. [7] The data do not support this evolution path: in contrast to 
the alignment models, the field decay model does not account for 
the increase in pulsewidth values for older pulsars. 

From this work we cannot set an upper limit on the magnetic 
field decay rate. Clearly, weak magnetic-field decay could still oc- 
cur concurrently with the alignment mechanism that causes the ob- 
served increase in pulsewidth for old pulsars. We can say that mag- 
netic field decay is not the dominant evolutionary factor: both the 
W50 and Wio fits are consistent with magnetic field decay not be- 
ing a significant factor in pulsar evolution. To model the combined 
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Table 3. Three-parameter least-squares fits of Models I, II and III. 



Rankin (1993) Data 





Tl -y 


Poo (deg) A^jog 




V 2 
*rcd 






ATNF Cat W 10 Fits (dof 


= 26) 




I 


2.47 


5.6 ±0.2 6.5 ±0.2 


0.71 ±0.08 


1.73 


JJ 




4.11 ±0.12 6.1 ±0.2 


0.83 ±0.08 


1 74 


TTT 
111 


Z.ZJ 


1.99 ±0.08 6.2 ±0.3 


0.89 ±0.08 


Y.lL 






ATNF Cat W 50 Fits (dof 


= 44) 








n 7 = 2.3 (no clear minimum in \ .) 




I 


2.3 


2.78 ±0.08 5.4 ±0.3 


1.02 ±0.08 


1.47 


II 


2.3 


2.10 ±0.07 4.8 ±0.4 


1.20 ±0.08 


1.46 


III 


2.3 


0.87 ±0.04 5.1 ±0.4 


1.26 ±0.09 


1.46 






Parkes MB W w Fits (dof 


= 10) 




I 


2.19 


5.5 ±0.3 6.1 ±0.4 


0.66 ±0.13 


1.05 


II 


2.30 


4.5 ±0.2 6.1 ±0.4 


0.71 ±0.14 


1.18 


III 


2.05 


1.84 ±0.15 5.6 ±0.5 


0.90 ±0.14 


1.16 






Parkes MB W 50 Fits (dof 


= 28) 




I 


2.48 


3.33 ±0.09 6.1 ±0.3 


0.83 ±0.09 


1.22 


II 


2.34 


2.36 ± 0.08 5.3 ±0.4 


1.08 ±0.09 


1.25 


III 


2.20 


0.87 ±0.05 5.1 ±0.4 


1.21 ±0.09 


1.23 
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Figure 7. The mean angular pulsewidths obtained from the magnetic field 
decay model, as fitted to binned ATNF Cat W\q data, are plotted against 
characteristic age. It has been assumed here that each pulsar has a constant 
inclination angle, oq. with the distribution of those angles determined by 
either Model I (dashed), II (solid) or III (dotted) of inclination angle as 
outlined in Section l3~3l the three curves are indistinguishable. 



effects of alignment and field decay is beyond the scope of this pa- 
per because of the relatively strong parameter cross-correlation that 
already exists in the alignment model alone. 



5.2 Inclination angle evolution 

To investigate the question of the a lignm ent of the magnetic 
an d spin axes. | Tauris & M anchester used the data sets 

of iRankinl jl993bh and iGouldl dl994l) that give estimates of a 
for a few hundred pulsars. These two datasets are derived using 
sli ghtly different methods and a ssumptions, as is well summarized 
bv lTauris & Manchester] |l998). Tauris & Manchester binned each 
data set in a to examine the (log t c ) (a) dependence. To these they 
fitted straight lines (their fig. 6, p. 632) in order to estimate the 
alignment time-scale, which they found to be around 10 7 yr. 
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Figure 8. Mean angle of inclination, (a), as a function of log charac- 
teristic age. The data are determined from models of beam morphology 
and polarization data and ar e taken from the two references considered by 
iTauris & Manchestejfl998T) . Binned data are shown with la error bars, and 
exclude low-field ( millisecond) pulsars having B SUI { ^ 4x 10 9 G. The top 
figure shows the Rankin ( 1993b, tables 2 and 4-8) inclination angles for 148 
pulsars; the central 3 bins each contain 30 p oints, while the first and last bins 
contain 29. The bottom figure shows the lGould jl994l tables 1-6) inclina- 
tion angles for 301 pulsars; of the 10 bins, all contain 30 points, except the 
sixth bin which has 31. The plotted curves are (o)(logt c ), equation i45\ , 
for each of the three models of initial alignment (I, II or III), using the pa- 
rameter values determined from the best fits to the ATNF Cat Wio data. 
The goodness-of-fit of the model curves to the binned data is measured us- 
ing the reduced chi-squared statistic, X^ ed , me va ' ues f° r eacn curve 
are shown on the plots. Models I and III are seen to give a poor fit to both 
the lRan kin and Goukj me an inclin ation-angle data, while Model II gives a 
particularly good fit to the Rankin data. 



Conversely, here the mean inclination angle associated with 
a log characteristic age is found by averaging the arcsin of equa- 
tion l !13t , for sin a, over the P(ao) and P(log r) distributions: 

<a) (log* e ) = 

oo it/2 

Iff,, . „ „ „ „, \ , dr 

A 



a(logt c )f b f t P(a )P(logr)da 



(45) 



Using the fitted-parameter values found in Section[4]from the vari- 
ous pulsewidth data sets, the resulting (a) (log t c ) curves for Mod- 
els I, II and III are plotted in Fig. [8] 

Also plotted are the data from iRankitj jl993bl) and iGould 
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Jl994h . in each case binned with about 30 points per bin. The error 
bars are lcr standard errors computed from the binning. The actual 
errors in each of the a values are difficult to evaluate, since they 
de pend on the validity of th e underlying assumptions, as discussed 
bv lManchester et all JT998). As with the pulsewidth data, low-field 
(millisecond) pulsars having _B sur f ^ 4 x 10 9 G have been removed 
from the data. 

The model curves are compared with the binned data using the 
reduced chi-squared statistic, x? e d> which is shown for each model 
curve in Fig. [8] It was found that the inclination angle data are in- 
consistent with Models I an d III, but compa tible with Model II. In- 
deed Model II predicts the Rankin ( 1993b) data particularly well, 
withaxLi ~ L 
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5.3 Beaming fraction evolution 

From equations d 1 3 1 ) and l |15l > for sin a(t c ) and p(t c ), equation d30t 
for ft becomes 



f b (t c ) =sina [l+4t c /(n- l)r]" 1/2 

x sin { Poo [1 + (n - l)r/4t c p /Cra - 1) } 



(46) 



The mean beaming fraction as a function of characteristic 
age, (/;,) (t c ), is a key quantity in pulsar population studies. It is 
found by integrating /(, (t c ) over the distribution functions P(«o) 
and P(log r). Using the parameter values from the best fit to the 
ATNF Cat Wio data set, (f b ) (t c ) is plotted in Fig.|9]for the three 
initial alignment models. 

Fig. |9] shows strong differences in the (/j) (t c ) among the 
three models for young pulsars because of the differing initial con- 
ditions, with gradual reduction to similar values for older pulsars. 
The figure implies that the smallness of the observed sample size 
for old pulsars is significantly influenced by a low probability of 
observation. Because of completeness issues, and luminosity evo- 
lution uncertainty, we cannot evaluate the relative significance of 
beaming fraction evolution and pulsar 'death' (through reducing 
radio emission) in determining the number of observable old pul- 
sars. 

For individual pulsars, the beaming fraction evolution will de- 
pend on the pulsar's ao and r 7 values, but the scaled beaming frac- 
tion fi, (t c )/ sin Qo is independent of ao- It is plotted in Fig. [9] for 
the Model II fit to the ATNF Cat Wio data, using various values of 
Tj. It is note-worthy that the scaled beaming fraction is consider- 
ably larger for pulsars with a longer alignment time-scale, and so 
these pulsars are more likely to be observed than those with a small 
alignment time-scale, even at young ages. 



5.4 Alignment time-scales 

From Table [3] the best-fitting mean log alignment time-scale, 
M-y.log, for Model II is found to be 6.1 ± 0.2 for the ATNF Cat 
Wio fit, 6.1 ± 0.4 for the Parkes MB W w fit, and 5.3 ± 0.4 for 
the Parkes MB W50 fit. The first and last of these three fits have 
the largest numbers of degrees of freedom, and they agree within 
two standard deviations of error. Furthermore, a spread in align- 
ment time-scales is needed to accurately model the data (CB86), 
and for the three fits just mentioned, the standard deviations, <7i og , 
are 0.83±0.08, 0.71±0.14 and 1.08±0.09 respectively, indicative 
of quite a large spread in alignment time-scales about the means. 

The mean time-scales quoted above are smaller than previous 
values obtained. The analyses of CB83, CB86 and C93 yielded a 
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Figure 9. (top) The mean beaming fraction, (fb), is plotted against log 
characteristic age for Models I, II and III, as labeled, (bottom) The scaled 
beaming fraction, / sin ao , is plotted against log characteristic age for 
Model II. The dashed curve has r 7 = /^ 7> i og , the dot-dashed curve has 
t 7 = £t 7i io g — o"j og , and the dotted curve has r 7 = At 7i i og + o"i og . For 
comparison, the corresponding Model II mean beaming fraction from the 
top plot is shown as a solid line. Both plots use the the ATNF Cat Wio 
best-fitting parameters, Table|3] 



Table 4. Alignment time-scales from different authors. 



Authors 



Alignment time-scale (yr) 



CB83, CB86, C93 

Lvne & Manchester (1988) 

Xu&Wu(1991) 

Kuz'min & Wu (1992) 

Tauris & Manchester (1998) 

Weltevrede & Johnston (2008) 



2 X 10 7 
10 7 
1.5 x 10 7 

2 x 10 7 
10 7 

7 x 10 7 



mean alignment time-scale of 2 x 10 7 yr. Other authors have ar- 
gued in favour of magnetic alignment using different methods and 
obtained similar al ignment time-scales of aro u nd 10 7 yr. See fo r 
examples Table | 4l l£vne & Manchesteil dl988h. Ixu & Wul Jl99lh . 
iKuz'min&Wul d 19921) andlTauris & Manchesteil 1 1998b all used 
polarization data, while IWeltevrede & Johnston! 1 20081) used inter- 
pulse pulsar statistics. 

From equation 1351 . the fitted value of /ii og is determined by 
the fitted value of ^i 7 ,i og and the value of 7: 



Miog = M7,iog - log (27) , (47) 
as plotted in Fig. [10] We see that a mean alignment time-scale 
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Figure 10. Plots of equation (47): the mean log alignment time-scale, /ii og , 
as a function of the parameter 7, which was introduced in equation (4)- The 
solid curve corresponds to £t 7l i g = 6.1 from the ATNF Cat Wio Model II 
fit (Table [3), and the dashed curves are three standard deviations of error 
away from the black curve, /U~ w = 6.1 ± 0.6. Previous studies have 
found a log alignment time-scale > 7. This is inconsistent here with the 
empirical finding that 7 1/2. 



10 w °b > 10 7 yr requires 7 > 1; that is inconsistent with previous 
findings that 7 is between 1/2 and 2/3. Conversely, for that range 
of 7-values, the curves in Fig.[l0]give a mean alignment time-scale 

: yr. 



5.5 Pulsar aging 

Caution is required in using the characteristic age as a measure of 
the true ages of pulsars, as will now be explored. Defining i 7 = 27t 
and using equations ( 122b and J24t involving r 7 , it follows from 
equation dl2t for t(t c ) that 



■ 111 



1 + 



4t e 



(jh 



It, 



(48) 



The mean log age, (logt 7 ), as a function of the log charac- 
teristic age, logt c , is found by integrating over the distribu- 
tion functions P(arj) and P(logr). Consider the best fit to the 
ATNF Cat Wio data having parameters as listed in Table[5] A plot 
of (log t 7 ) (logf c ) for Model II is shown in Fig. ll H it shows that 
that (logi 7 ) and logic agree up to logi c ~ 10^ 7l i og , after which 
time (logfcy) starts to fall below logic. The curves for Models I 
and III are similar. 

However, the mean value (log t) (logic) applies to the popu- 
lation as a whole. For individual pulsars, the plots of logt(logt c ) 
depend on the value of t 7 for that pulsar. Various tracks for in- 
dividual pulsars are shown in Fig.[TT] it is evident that in general 
the characteristic age diverges much more rapidly from the true age 
than for the populations mean ages, with the divergence being more 
rapid the smaller the alignment time-scale. Recall, however, that 
pulsars with shorter alignment time-scales have a smaller beaming 
fraction than those with longer time-scales, and so are less likely to 
be observed. 

This highlights the uncertainty that may be present in previous 
studies that have used the characteristic age as a measure of true 
age when attempting to measure the alignment time-scales. The 
characteristic age can over-estimate the true age, and so will tend 
to give larger alignment time-scales if it is the time measure. 
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Figure 11. Plots of mean log age, (logi 7 ), (solid curve) and log age 


logi 7 


(broken curves) as 


functions of log characteristic 


age for Model II 



(dashed) r 7 = /i 7 i og , (dot-dashed) r 7 = £t 7i i g — ""log^ an d (dotted) 
r 7 = M-y^og + ci g. The parameters are from the best fit to the ATNF Cat 
Wio data set (Table[5J, and the straight line is log t 7 = log t c . 



6 CONCLUSIONS 

From Fig.(6] showing mean pulsewidths versus log t c , we have seen 
that the pulsewidth data clearly favour pulsar evolution through 
magnetic alignment, and not magnetic field decay alone nor pro- 
gressive counter-alignment. There is a clear minimum of mean 
pulsewidth at about t c w 10 6 ' 9 yr. The alignment process causes 
a significant increase in the mean Wio pulsewidth for all charac- 
teristic age bins greater than 3 x 10 7 yr. This is particularly signifi- 
cant in view of the fact that the ATNF catalogue contains numerous 
pulsars with t c > 10 s yr, which were not available to the earlier 
studies. 

The three models of initial alignment that were specified in 
Section [33l cannot be tested by fitting the Candy-Blair models to 
the pulsewidth data alone. However, once a Candy-Blair model has 
been fitted to the pulsewidth data, it can be used to predict mean 
inclination angle evolu tion, (a) (log i e ) ■ Co mparing these predic- 
tions to the datasets of Rankin ( 1993b) and lGould jl994l) . it was 
found that only Model II (a random initial inclination angle) gives 
a good fit to the data, as shown in Fig. [8] 

The alignment time-scale found from fitting Model II to the 
ATNF Cat Wio data (see Tables □ and [3} is 10 6 1 yr, which is at 
least a factor of ten shorter than the timescales determined by a 
number of other authors (see Tabled- The fit requires that the log 
alignment time-scale varies with a standard deviation of about 0.83 
across the observed population. The fit gives a limiting cone open- 
ing half-angle of 4? 11, and a braking index parameter n 7 of 2.30, 
consistent with the mean, 2.42 ± 0.05, of the measured apparent 
braking indices (equations (O and d25t , and Table O. Consistent 
with empirical observations, the fit assumes the parameter 7 = 1/2 
(see equation (Q). For other values of 7 the fitted parameters n and 
/xio g are determined from equations j_ 1 b and {35} respectively. A 
plot of the mean beaming fraction for this fit is given in Fig. [9] 

The initial predictions of pulsar alignment due to electromag- 
netic braking jDavis & Goldsteinlll97dlMichel & Goldwirelll970l) 
predicted time-scales of around 10 3 -10 4 yr, which was rather 
small. The time-scales will be increased if pulsars are born with 
their magnetic and rotation axes nearly orthogonal. Ijonesl J 19761 : 
Il976al ; Il976bl) found that the inclusion of a decreasing dissipative 
torque initially brought all pulsars close to counter-alignment, from 
which time the normal electromagnetic decay of a would begin. 



12 M.D.T. Young, L. S. Chan, R. R. Burman and D. G. Blair 



Jones predicted that the alignment phase started at around 10 3 - 
10 4 yr, with a time-scale of around 10 6 yr, which was believed to 
be more in accord with the observed population. 

We have found alignment time-scales consistent with the pre- 
dictions of ljonesl though with a significant spread in values, down 
to the time-scales consistent with the absence of the counter- 
alignment phase. Contrary to the model of Jones, though, we do 
not find evidence for all pulsars starting their alignment phase with 
Qo close to 90° (i.e. Model I). This interesting counter-point to the 
model of Jones deserves further investigation. 
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APPENDIX A: PULSEWIDTH AVERAGING INTEGRAL 

The integral that averages the pulsewidth over a random distribu- 
tion of viewing angles can be evaluated analytically, as follows. As 
above, a is the pulsar magnetic-spin inclination angle, p the emis- 
sion cone half-width, and ( is the angle between the observer's di- 
rection and the pulsar's spin axis. 



Al Observer dependence of pulsewidth 

For an individual pulsar, the variation of its observed angular 
pulsewidth with observer v iewing angle is given by the formula 
ilManchester & Tavloill 19771. p. 218) 



W(a,p, = 2arccos 



cos p — cos a cos £ 
sin a sin£ 



(Al) 



which is equation i26l in the main text above. 

In considering the variation of W with viewing angle, it is 
helpful to employ the variable x = cos (. Equation dAlt takes the 
form 



W (a, p,x) = 2 arccos[/(a;)] 

with the abbreviations 

1 b — ax 



a = cos a , and b = cos p . 



(A2) 



(A3) 



sin a (1 — x 2 ) 1 / 2 
Note that 

(sin a) /'(a;) = (bx - a) /{I - x 2 f /2 and (A4) 
(sin a) f"(x) = (6 - Sax + 2bx 2 )/(l - x 2 f' 2 . (A5) 

The function f(x), which is just the cosine of the pulse angu- 
lar half-width, is < 1 throughout the cut across the emission beam 
by the sight line to any 'illuminated' observer. It reaches 1 in mag- 
nitude (i.e. f 2 = 1) at the pulse visibility limits £ = | a =p p\, 
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corresponding to x = x±, where 

x± = cos(q =p p) — ab ± A 
A = (1 -a 2 ) 1/2 (l -b 2 ) 1 ' 2 



with 



Substituting dA6t into dA3t gives 



/(*±) 



6(1 



2x1/2 



To(l 



[1- (ab± A) 2 ] 1 /2 



(A6) 
(A7) 



(A8) 



This shows that at a;_, the 'outer' visibility limit, / is always posi- 
tive, i.e. +1. At x+, the 'inner' visibility limit, / is +1 for a > p, 
for a = p, and — 1 for p > a. Let's visualize these three cases in 
turn. 

(1) For a > p, the visibility limits x± both correspond to zero 
pulsewidth, W(x) = 0. Equations dA4t and dA5t show that the 
minimum value (b 2 — a 2 ) 1 ' 2 /(sin a) of f(x) - corresponding to 
the greatest pulsewidth - occurs for the cone of observers at x = 

a/b: 



2 arccos 



(sinp) 2 /(sina) 2 ] 1/2 | 



(A9) 



at cos ( = (cos a)/(cos p). 

(2) For a = p, one can visualize the emission cone as rolling 
around the spin axis, which is then the inner visibility limit. In this 
case, f{x) = (cota)[(l - x)/(l + x)] 1/2 ; also x+ = +1 (ob- 
server looking along the spin axis) so f(x+) = 0, corresponding 
to W(x+) = tv. So W(x) increases monotonically from at the 
outer visibility limit to tt on the spin axis in this special case. 

(3) For p > a, the beam will be 'always on' (unpulsed) to ob- 
servers within (p — a) of the spin axis. So the pulse visibility limit 
x+ then corresponds to W — 2n, with f(x) ranging from +1 at 
X- to — 1 at x+, with no minimum: W(x) increases monotonically 
from at the outer visibility limit to 27r at the inner one. 



A2 Evaluating the averaging integral 

Integrating equation dAlt for W(a,p, over a random distribu- 
tion of viewing angles, using the observer distribution function 
(sin C)/(2/f>), gives the mean angular pulsewidth, for pulsars of 
given a and p: 

(W)(a,p) = (2ft)- 1 [ V(p,a,C)sinCdC (A10) 
= /(,~ / arccos[/(a;)] da; . (All) 



It is integration of the observer distribution function (sin £) /2 over 
the same range that yields the beaming factor ft = sin a sin p. 
Note for use below that 

(sin a) [l-/ 2 (:r)] 1/2 = G 1/2 (x)/(l - a; 2 ) 1/2 , (A12) 
With G(x) = (x-x-)(x+-x) . (A13) 

The function G(x) describes an inverted parabola: it is zero at the 
limits of integration x± and positive in between, with a maximum 
value (1 — a 2 )(l — b 2 ) = sin 2 a sin 2 p (= ft' 2 ) at x = ab = 
cos a cos p. 

Applying integration by parts in equation \Al It for (W) re- 
duces the problem to one of evaluating an integral of an algebraic 
function. Using d(arccosX) = —(1 — X 2 )~ 1 ^ 2 dX together with 
equations JA4b and JA12b for df/dx and (1 — f 2 ) 1 ^ 2 gives 



+ 



+ x(bx — a) dx 
(l-x 2 )G 1 / 2 (x) 



(A14) 



because the integrated part x arccos [/(a;)] - which by equa- 
tion iA2l is just x W(x) /2- always vanishes at the limit X-, where 
/ = +1 always. The first term on the RHS of equation dA14t is 

for a>p \f(x+) = +l], 
ty/2 for a = p [f(x+) = , x+ = +1} , (A15) 
7rcos(p — a) for a < p [/(a;+) = —1] . 

Re-arranging the integrand in equation dA14t gives 

(W)(a,p)f b =^±W(x + ) 



+ 



b — a b + a \ dx 

+ — r — o 



2(1 -x) 2(1 + ac) 



G 1 / 2 (a 



(A 16) 



This result shows that we can employ the following indefinite 
integrals, in which A and G are defined by equations dA7t and 
dAL3l : 



G~ 1/2 (x)dx = , 



da; 



i[(a; - ab)/A] and 



da; 



(A17) 



(l±x)G 1 ^(x) J (l± x )[(x-x.)(x+-x)] 1/2 
2 



■ arctan 



(x — a;_)(l ± x+) 



\a±b\ V (x+ — x)(l ± X-) 

where x± = ab ± A and equation dA7t for A have been used 



(A18) 



So the required definite integrals are 

G~ 1 ^ 2 (a;) da; = rr , and 
da; 7r 



(l±x)G 1 / 2 (x) 



\a±b\ 



(A19) 



(A20) 



Thanks to these definite integrals, equation dA 1 6b for the observer- 
averaged pulsewidth evaluates to 

(W)(a,p)f b = ^±W(x+) + ^[ S ign(b-a) + l-2b] . (A21) 
That is: 

{W){a,p)f b = 

[l-cosp]7r for a^p, 
[cos(p — a) — cosp] 7r for a^f, 

as used in the main text - equation d38t above. 



(W)(a,p)f b = x+ arccos[/(aj+)] 



